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Adiabatic processes in the quantum Ising model and the anisotropic Heisenberg model are dis- 
cussed. The adiabatic processes are assumed to consist in the slow variation of the strength of 
the magnetic field that environs the spin-systems. These processes are of current interest in the 
treatment of cold atoms in optical lattices and in Adiabatic Quantum Computation. We determine 
the probability that, during an adiabatic passage starting from the ground state, states with higher 
energy are excited. 
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I. INTRODUCTION 

The analysis of ground state properties of spin sys- 
tems in the vicinity of quantum phase transitions is a 
contemporary issue in statistical physics and condensed 
matter physics. These systems display a rich variety of 
phenomena which explain many exotic properties of ma- 
terials at very low temperature. Very recently, it has 
been recognized that cold atoms in optical lattices can 
be very well described in terms of certain spin Hamilto- 
nians These systems can be very well controlled at 

the quantum level, which may allow us to observe physi- 
cal phenomena predicted for spin systems and that have, 
so far, not been observed with other systems. Moreover, 
they may allow to build a quantum simulator which may 
shed some light in unresolved issues in condensed matter 
physics 13. 

In the strong interacting regime and at sufficiently low 
temperature, atoms confined in optical lattices tend to 
distribute in the so-called Mott phase; that is, under the 
appropriate conditions each lattice site is occupied by 
a single atom. Each atom still interacts with the near- 
est neighbor via virtual tunneling, which gives rise to an 
effective interaction between the internal levels of neigh- 
boring atoms. Under certain conditions, these interac- 
tions can be described in terms of the Ising Model or the 
anisotropic Heisenberg Model with a transverse magnetic 
field 1^, ^] . Changing the laser parameters amounts to 
varying the parameters of these models. Both models 
display quantum phase transitions '7| for certain values 
of these parameters. Thus, with these systems it is possi- 
ble to investigate the dynamics of quantum phase transi- 
tions by, for example, adiabatically changing these laser 
properties |24j . An important question in this case is to 
what extent the process can remain adiabatic since near 
a phase transtion the energy levels tend to get closer as 
the number of particles increases |3- In this paper we 
will study this problem for both Hamiltonians. For the 
exactly solvable Ising Model, the adiabatic process will 
be investigated in detail and analytical results for the ex- 
citation probability will be given. For the yet unsolved 
anisotropic Heisenberg Model, the adiabatic process will 
be investigated by means of perturbation theory. 

Another subject in which the results about adiabatic 



processes in spin-systems are relevant is Adiabatic Quan- 
tum Computation |^ Il0j| . In Adiabatic Quantum 
Computation, solutions to mathematical and physical 
problems are obtained by simulating adiabatic processes 
on quantum computers [III IT^ ITsl IT^ . An interest- 
ing problem related to Adiabatic Quantum Computation 
is the investigation of the ground state of spin-systems. 
This problem is known to be difficult to be solved on clas- 
sical computers, since the required resources in time and 
space scale exponentially with the number of particles the 
system consists of. Because of this exponential scaling, 
the ground state properties can only be obtained for spin- 
systems consisting of a few dozens of particles. In addi- 
tion to that, analytical calculations are rarely possible, 
such that ground state properties of many spin-systems 
are still unknown. The question whether the ground 
state of spin-systems can be investigated efficiently by 
means of Adiabatic Quantum Computation is related to 
the question how slowly parameters of the system must 
be changed such that processes remain adiabatic. This 
question will be answered in detail for the quantum Ising 
model. Even though this model is exactly solvable and its 
ground state is known, the results about this model shed 
light on the efficiency of adiabatic quantum algorithms 
investigating the ground state of more complicated spin- 
systems. In the case of the anisotropic Heisenberg model, 
statements about the efficiency will be made in regimes 
amenable to perturbative treatment. 

The investigation of the ground state of spin-systems 
by Adiabatic Quantum Computation is an example for 
a quantum algorithm that is feasible with current tech- 
nology: This algorithm only requires a low number of 
qubits, of the order of 50, in order to exceed classical 
computations. In addition, as we will show, the algo- 
rithm is sufficiently robust against errors and imperfec- 
tions, such that no quantum error correction codes are 
necessary. On top of that, the algorithm can be imple- 
mented by any experimental realization of a Universal 
Quantum Simulator 12], like that based on optical lat- 
tices and arrays of microtraps 3]. 

The paper is outlined as follows: In the first part, gen- 
eral statements about the quantum Ising model and the 
anisotropic Heisenberg model are made. In the second 
part, the investigation of the ground state of spin-systems 
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by Adiabatic Quantum Computation is discussed. In the 
third and forth part, adiabatic processes in the quantum 
Ising model and the anisotropic Heisenberg model are 
studied. The effects of perturbations on these adiabatic 
processes are discussed in the last part. 



II. PRELIMINARIES 



A. Ising Model 

The quantum Ising chain 0, ^3 a perfectly suited 
model for investigating adiabatic processes because it is 
exactly solvable and shows a quantum phase transition. 
The Hamiltonian that describes this model reads 



H 



N-l 



i=0 



(1) 



I.e. 



Gy^ and a^z^ describe the spin of 



The chain is assumed to be cyclic 
The Pauli operators 

the I*'' particle in the chain. The dimensionless variable g 
determines the strength of the transverse magnetic field. 
The parameter J is considered as a positive quantity that 
fixes the microscopic energy scale. N denotes the number 
of particles in the chain. 

The diagonalization of Hamiltonian J^l can be 
formed by means of Jordan- Wigner transformation 

This transformation maps the Pauli operators (Tx , uy 

and af^ on Fermi operators. The details of the diago- 
nalization procedure can be gathered from appendix 1X1 
After diagonalization, the Hamiltonian of the quantum 
Ising model assumes a simple form, namely, 



3er- 



H/J = gN 



A. 



N even 
N odd 



with N denoting the number of fermions. The variables 
Ak and Ak denote energies of single fermions and the 77's 
are fermionic operators. 

The eigenvalues of Hamiltonian Q as functions of the 
transverse field g are visualized in figure ^ It is eas- 
ily shown that the ground state energy as a function 
of g is non-analytic at the point g = 1 in the thermo- 
dynamic limit. This non-analytic point is an indication 
for a quantum phase transition at the quantum crit- 
ical point g — 1. This quantum phase transition sepa- 
rates two phases at zero temperature: an ordered phase 
for 5 < 1 and a disordered phase for g > 1. In the or- 
dered phase, the interaction between the particles aligns 
the spins parallel or antiparallel to the crystalline axes. 
In the disordered phase, the external magnetic field in- 
duces spin-flips that destroy this order. 

The symmetries that Hamiltonian possesses are the 
Z2-symmetry and the translational symmetry. The Z2- 
symmetry corresponds to a reflection of the spin-vectors 
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FIG. 1: Spectrum of an Ising chain consisting of 5 particles 
{N = 5) as a function of the transverse magnetic field g. The 
dotted line marks the critical point at which the quantum 
phase transition takes place. 



at the x-axes. It is generated by the operator 



N-l 

n- 

1=0 



This operator possesses two different eigenvalues, 
1 and —1, such that eigenstates can be classified by 
two different Z2-symmetries. In terms of Fermi opera- 
tors, states with an even number of fermions possess Z2- 
symmetry 1 and states with an odd number of fermions 
possess Z2-symmetry —1. The translational symmetry 
is generated by the translation-operator which has the 
property to right-shift all product states: 

T\zi) ■ ■ - {zn) ^ \ zn)\zo) ■ ■ ■{ ZN^i ) 

This operator possesses N different eigenvalues e"*^'^"/^ 
(n = 0, . . . , A^ — 1), such that the eigenstates can be clas- 
sified by A^ different translational symmetries. 



B. Heisenberg Model 

The anisotropic Heisenberg chain [T7LI1MIT9I I is charac- 
terized by anisotropic internal interactions. The Hamil- 
tonian that defines this model reads 

H^-JJ2 {A^cri^V^D + A,a«a('+i) + (2) 

The chain is assumed to be cyclic. The properties of 
the interactions between the particles of the chain are 
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FIG. 2: Spectrum of the anisotropic Heisenberg model as 
a function of the external magnetic field g. Assumptions: 
iV = 5, = 0.1, Aj, = 0.3 and A^ = 1; The table on the 
right hand side gives informations about the degeneracy D, 
the group-index n and the translational-symmetry k of each 
level. 



described by the anisotropy-parameters A^;, Aj, and A^ 
and the strength of the external magnetic field is deter- 
mined by the parameter g. 

The anisotropic Heisenberg model possesses the same 
symmetries as the quantum Ising model: the Z2- 
symmetry and the translational symmetry. Thus, eigen- 
states can be classified by two different Z2-symmetries 
(1 and —1) and by N different translational symmetries 
(e-'^^"/^ ^ith n = 0, . . . , TV - 1). 

Another possibility to classify the eigenstates is to clas- 
sify them according to their behavior in a very strong 
magnetic field. In figure[21the energies of the eigenstates 
are plotted as functions of the field-strength g. From 
this figure it can be gathered that the eigenstates clus- 
ter in -I- 1 groups because of the Zeeman-shifts. All 
states in each group have the property that they tend to 
the same state in the limit g — > 00. In the n^^ group, 
for example, all states tend to a state that has all but 
n spins aligned in x-direction and n spins pointing in 
— a;-direction. All eigenstates in the n}^ group can be 
classified by a group-index n. This group-index is re- 
lated to the Z2-symmetry of the states. It can be shown 
that states with an even group-index possess the Z2- 
symmetry 1 and states with an odd group-index possess 
the Z2-symmetry —1. 

In order to take into account that several eigenstates 
with group-index n and translational symmetry k may 
exist, a third index d is required to distinguish eigenstates 
with equal values for n and k. In summary, all eigenstates 
of H can be identified by the three indices n, k and d. 



They will be denoted as 

I i'nkd ) , 

where n e {0, . . . , TV} and fc e {0, . . . , A^ - 1}. Evidently, 

the kets | V'nfcd ) are common eigenstates of H , Z2 and T . 

III. QUANTUM SIMULATIONS OF SPIN 
SYSTEMS 

Before studying adiabatic processes in spin-systems, a 
relevant application shall be discussed: the investigation 
of the ground state of spin-systems with a quantum com- 
puter. This investigation could be performed by means 
of Adiabatic Quantum Computation 0, 0, 0| . 

The scheme of Adiabatic Quantum Computation reads 
as follows: First, the quantum computer is prepared in 
the ground state of a simple beginning Hamiltonian, the 
ground state of which is known. Then, using the Univer- 
sal Quantum Simulator 12], a time evolution according 
to a time-dependent Hamiltonian H{t) is simulated which 
adiabatically interpolates between the simple beginning 
Hamiltonian and a complicated problem Hamiltonian, 
the ground state of which shall be investigated. Because 
of the adiabatic variation, the quantum computer always 
stays in the ground state of the time-dependent Hamil- 
tonian H{t) and is finally prepared in the ground state 
of the complicated problem Hamiltonian. Finally, by ap- 
propriate measurements, information about the ground 
state of the problem Hamiltonian is obtained. 

Evidently, in the case of spin-systems, the problem 
Hamiltonian is the Hamiltonian that describes the spin- 
system of interest. The ground state properties of this 
Hamiltonian are usually interesting as functions of the 
strength g of an external magnetic field. Thus, it is obvi- 
ous to lay out the algorithm such that the field-strength g 
is adiabatically varied and tunes between the simple be- 
ginning Hamiltonian and the problem Hamiltonian. The 
beginning Hamiltonian can either be the Hamiltonian of 
the spin-system for g — or the Hamiltonian of the spin- 
system for (7 — > 00. The algorithm then consists either 
in preparing the quantum computer in the ground state 
for g = and adiabatically increasing the field-strength g 
or in preparing the quantum computer in the ground 
state for g —^ 00 and adiabatically decreasing the field- 
strength. The ground state for g — > cxd is simple since it 
has all spins aligned in the direction of the external field. 
Whether the ground state for g — can be prepared or 
not depends on the spin-system under study. In fact, 
many models exist that are solvable for g = and, in the 
case of these models, g = can be chosen as a starting 
point. 

The duration T of the algorithm is related to the 
change rate of the parameter g. The parameter g must 
be changed slowly enough, such that the time evolu- 
tion remains adiabatic. The change rate can be deter- 
mined mathematically by means of the Adiabatic The- 
orem [23, 1^ [22]. This theorem deals with the solu- 
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tion of the Schrodinger equation in the case of a time- 
dependent and slowly varying Hamiltonian. The state- 
ments of the Adiabatic Theorem are the following: In 
the limit T — > cxd, the system always stays in the ground 
state. Thus, eigenstates of the beginning Hamiltonian 
are mapped with certainty on eigenstates of the problem 
Hamiltonian and there is only a change in the phase. In 
reality, however, the duration T of the time evolution 
is finite. Thus, the question that has to be answered 
is how large the duration T must be, such that, with a 
high probability, the system stays in the ground state. In 
other words, the probability that eigenstates with higher 
energy are excited must be negligible. 

A rough criterion on the duration T that guarantees 
that the excitation probability is negligible reads 

r»A. (3) 

£ is thereby a quantity that depends on the derivative of 
the Hamiltonian H{s) with respect to s and A denotes 
the minimum energy difference between the ground state 
and the first excited state of H{s). Usually, £ scales 
polynomially with the number of particles, such that the 
efficiency mainly depends on the behavior of the mini- 
mum energy difference A as a function of the particle 
number N. The energy difference usually reaches its 
minimum at an avoided level-crossing. At an avoided 
level-crossing, the ground state and the first excited state 
approach as the number of particles increases, such that 
the duration of the algorithm will always increase as the 
number of particles increases. The efficiency of the algo- 
rithm is then determined by the velocity with that the 
first two states approach. If, on the one hand, the first 
two states approach exponentially fast with a growing 
particle number, the algorithm is not efficient. If, on the 
other hand, the two states approach polynomially fast, 
the algorithm can be considered to be efficient. 

What is seen from this discussion is that knowledge of 
the spectrum of the Hamiltonian H(s) is required in order 
to make statements about the duration of the algorithm. 
However, the spectrum is not known, in general. If the 
spectrum was known, the ground state of the problem 
Hamiltonian would be known, as well, and no quantum 
computer would be required to calculate it. Thus, the 
only thing that can be done is to test the quantum algo- 
rithm on spin-systems that are solvable and see whether 
it is efficient or not, or to try to find estimations of the 
spectrum and use them to approximate the duration of 
the algorithm. 

Because of these difficulties in the mathematical de- 
termination of the duration of the algorithm, a simple 
experimental method is desirable that can be used to 
check whether a chosen change-rate is sufficient or not. 
An experimental method could look like this: First, the 
quantum computer is prepared in the ground state of the 
beginning Hamiltonian, e.g. the Hamiltonian of the spin- 
system for g = 0. Then, the field-strength g is increased 
with a chosen change rate up to a desired value and de- 
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FIG. 3: Plot of the excitation energies of the states | n ) as 
functions of g. The assumed number of particles A'^ is 51. 

creased again with the same change rate. At the end, a 
measurement in the Eigenbasis of the beginning Hamilto- 
nian is performed (which is known). From this measure- 
ment it can be deduced whether the quantum computer 
is still in the ground state or whether levels with higher 
energy have been excited. If the quantum computer is 
still in the ground state, the time evolution was adia- 
batic and the change rate was low enough. Otherwise, 
the change rate must be decreased and the experimental 
check must be performed once again. 

IV. ADIABATIC PROCESS IN THE ISING 
MODEL 

The adiabatic process under study consists in the slow 
variation of the strength of the magnetic field that envi- 
rons the Ising chain, i.e. the slow variation of parameter g 
in Hamiltonian (^) . This variation leads to excitations of 
states that possess the same symmetry as the initial state 
of the system (which is assumed to be the ground state) . 

The symmetries of the problem are, as discussed in 
section lll Al the Z2-symmetry and the translational sym- 
metry. Z2-symmetry imposes that only states consisting 
of an even number of fermions are excited. Translational 
symmetry makes further restrictions. The states that are 
excited in first order can be identified by means of the 
Adiabatic Approximation |2^. They read 

where n — 1, . . . , [^^^^J. The energies of these states 
with respect to the ground state energy are plotted in 
figure 13 as functions of the field-strength g. This figure 
shows that the energy difference between the lowest ex- 
cited state and the ground state tends to a minimum at 
the critical point g — 1- The minimum energy difference 
thereby amounts to ^ for N ^ 1. This observation al- 
ready allows to make a qualitative statement about the 
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FIG. 4: Excitation probability po-»n of the levels | n ) in 
the case of 501 particles (A*" = 501). The magnetic field is 
assumed to be varied from go = 5 to gi with change rate 
g equal to —0.0001. The indices n are reshufHed such that 
states with lower index have lower energy. 



duration of an adiabatic time evolution in the vicinity of 
the critical point: Since, according to formula (|3J), the 
duration scales with where A is the minimal energy 
difference between the lowest exited state and the ground 
state, the duration can be expected to scale with N^, i.e. 
with the square of the number of particles. This has the 
consequence that, for example, the investigation of the 
ground state of the quantum Ising model in the vicin- 
ity of the critical point by means of Adiabatic Quantum 
Computation will be efficient. 

These qualitative results can be clarified by calculat- 
ing the excitation probabilities pq — 

for all states In). 

This can be done analytically by making use of the Adi- 
abatic Approximation |2^. The details of the analytical 
calculation can be gathered from appendix IbI Another 
way to calculate the excitation probabilities po-»n con- 
sists in numerically solving the Heisenberg equations for 
the Fermi-operators rjk- Using this solution, it is simple 
to calculate expectation values of arbitrary observables 
after the adiabatic process. If it is assumed that the fi- 
nal state of the system is a superposition of the ground 
state and the 2-fermion states | n ) only, the excitation 
probabilities po->n can be calculated as well, since they 
are expressible in terms of expectation values of the op- 
erators 77lj77^„7y_„77„. 

In figure 01 the excitation probabilities po->n are plot- 
ted for an Ising chain consisting of 501 particles. The 
adiabatic process is assumed to consist in the linear de- 
crease of the magnetic field-strength g with change rate 
g = -0.0001 from go = 5 to 51 = 1.08, 1.05, 1.03 and 



FIG. 5: Excitation probability pE as a function of the final 
field strength gi. Assumptions: 51 particles (A'^ = 51); linear 
variation of g from go = 5 to gi with change rate g = —0.01. 
Solid line: Sum of all terms po^n', Dotted line: Numerical 
results (see text); Dashed lines: Analytical estimations (for- 
mulas @ and JSJl). 



1.0 respectively. It can be observed that, if gi is in the 
vicinity of the critical point, the lowest state is excited 
with maximal probability. 

The overall excitation probability is obtained by 
summing the excitation probabilities pq — of the respec- 
tive states. Analytical estimations for the overall excita- 
tion probability are obtained in three different regimes: 

(1) go » 1, 51 > 1 and N{gi - 1) > 1 

(2) go > 1, 51 > 1 and N{gi - 1) < 1 

(3) go > 1 and gi < 1 

• In regime (1), the analytical expression for reads 



Pe 



< 



{90-91? 



28 J2 



T2 [gi-iy 



■iV, 



(4) 



such that the duration T of the adiabatic process 
will scale with the square root of the number of 
particles N . 

• In regime (2), the overall excitation probability pe 
can be written as 

ft' (30 - 51)' 



PE 



< 



267^4 J2 T2 



-(4-3gi)A^\ 



(5) 



Thus, the duration T of the adiabatic process will 
show the previously mentioned A^^-scaling. 

• In regime (3), the overall excitation probability 
is constant as a function of gi (for a fixed change- 
rate g EE (go - gi)/T), namely, 

ft' (50 - 91? 



PE 



< 



267^4 J2 T2 



(6) 



and the duration T of the adiabatic process will 
show a quadratic scaling with the number of parti- 
cles N, just like in regime (2). 
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FIG. 6: Duration T of the adiabatic process as a function of 
the particle- number A'^. Assumptions: g is Unearly decreased 
from go = 5 to (?i = 0; the change-range g is adjusted such 
that the excitation probabihty pE always amounts to 0.05; 
Grosses: numerical results (see text); Dash-dotted line: an- 
alytical result i.e. T ft; 0.2832A'^'^. Diamonds: results 
presented in for the free-end chain. 



FIG. 7: Expectation value of H relative to the width of the 
spectrum, i.e. {{H) — Eo)/AE, as a function of the change- 
rate g. Assumptions: N = 51, go = 5 and gi — 0; solid 
line: numerical results (see text); dash-dotted line: indicates 
variance of H; dots: results obtained using the Adiabatic 
Approximation. 



In figure the overall excitation probability pE after 
the decrease of g from go = 2 down to gi is shown as a 
function of the final field-strength gi . The change rate g 
is assumed to be fixed at —0.01 and the number of parti- 
cles is taken to be 51. The solid line represents the sum 
of the terms po—>n that were calculated by means the 
Adiabatic Approximation. The dashed lines represent 
the analytical estimations I0J and ||SJ|. The dots repre- 
sent numerical results obtained by solving the Heisenberg 
equations for the Fermi operators. What can be observed 
is that the excitation probability increases fundamentally 
as the critical point g = 1 is crossed. This behavior is 
due to the fundamental approach of the ground state and 
the lowest excited state at the critical point. 

The duration T as a function of the number of par- 
ticles N of an adiabatic process consisting in the linear 
decrease of the field-strength g from go = 5 to gi = 
can be gathered from figure IHI The change-rate g is 
thereby adjusted such that the overall excitation prob- 
ability pe is always equal to 0.05. The analytic result © 
is represented by the dash-dotted line and the results ob- 
tained from numerically solving the Heisenberg equations 
for the ?7fe's are represented by crosses. In both cases, a 
quadratic scaling of the duration T with the number of 
particles N can be observed. 

Up to now it has been assumed that the field-strength 
is changed sufficiently slowly, such that the Adiabatic Ap- 
proximation is valid. Outside the adiabatic regime, 
states with a fermion-number 4, 6, 8, . . . are excited, ad- 
ditionally to the 2-fermion states | n ) . These states con- 
tribute to the overall excitation probability and must be 



taken into account. However, even outside the adiabatic 
regime, the mean energy of the system still remains very 
low. This can be gathered from figure [3 In this figure, 
the mean energy after the decrease of the field-strength g 
from go = 5 down to gi = is plotted as a function of the 
change-rate g for a 51-particle chain. The results about 
the mean energy were obtained by numerically solving 
the Heisenberg equations for the rj^s. As a consequence, 
processes marginally outside the adiabatic regime pre- 
pare the system in a state with very low temperature the 
properties of which are still interesting to be investigated 
with a quantum computer. 



V. ADIABATIC PROCESS IN THE 
HEISENBERG MODEL 

The adiabatic process in the anisotropic Heisenberg 
model is treated similarly as in the previous section: 
Symmetries are used to determine the states that are 
excited during the adiabatic process. Since the adiabatic 
process starts from the ground state, the excited states 
can be identified as the states that possess the same sym- 
metries as the ground state. Thus, using the notation of 
section IIIBI the excited states are the states with an 
even group-index n and translational symmetry k — 0. 
In practice, only states from the 2'"* group are excited in 
first order, namely, 

I '4'20d ) 

with d ranging from 1 to -^^j"'" is assumed to be odd). 
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FIG. 8: Excitation probability pe as a function of the change 
rate a = (f/i — go)/T in the case of 5- and 9-particle chains. 
Assumptions: /S.x = 0.1, Ay = 0.3, = 1; variation of g 
from go — 10 to gi — 5; Dotted line: Excitation probability 
obtained from a numerical simulation. Solid line: Approxi- 
mation (|7| of the excitation probability in the limit A'^ S> 1. 
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FIG. 9: Excitation probability pe as a function of the num- 
ber of particles N in the case of change rates a = —6 and 
a = —10. Assumptions: A^, — 0.1, A^ — 0.3, A^ = 1; vari- 
ation of g from go = 10 to gi = 5; Dotted line: Excitation 
probability obtain from a numerical simulation. Solid line: 
Approximation O of the excitation probability in the limit 
A > 1. 



Using the previous considerations, it is straightforward 
to study the efficiency of the algorithm in the case of the 
anisotropic Heisenberg model. What is required to be 
done is to find perturbative expansions of the ground 
state and the possible excited states and to use the Adi- 
abatic Approximation 23J to calculate the excitation 
probabilities. Of course, in order to be allowed to trun- 
cate the perturbative expansions after a few lower-order 
terms, it must be demand that the Hamiltonian consists 
of parts that are of a different order of magnitude. Be- 
cause of this, it is assumed that the strength g of the 
external magnetic field is always much larger than the 
strength of the internal interactions. The adiabatic pro- 
cess therefore consists in adiabatically decreasing the field 
strength g from go oo down to a certain value gi , which 
is still much larger than the strength of the internal in- 
teractions {gi ^ A^). 

The results that are obtained in this way are very simi- 
lar to results obtained for the quantum Ising model. De- 
tails about the calculations can be gathered from ap- 
pendix [CI The analytical expression that is obtained for 
the overall excitation probability pE reads 

PE<^iA,~Ay) — -giV, (7) 

such that the duration T of the adiabatic process will 



scale with the square root of the number of particles N. 
As a consequence, it will be possible to investigate the 
ground state of the anisotropic Heisenberg model effi- 
ciently by means of Adiabatic Quantum Computation 
in the regime where the external magnetic field is much 
stronger than the internal interactions. 

The analytical expression Q can be compared to the 
exact result (@J) presented in section Hvl for the quantum 
Ising model. As it is easily verified, the exact result is 
equal to expression Q) in the limit 3> 1 if the param- 
eters Aa;, Ay and are specified to A^; = 0, Ay = 
and Az = 1. 

Expression (0 can also be compared to results ob- 
tained from numerical simulations. This is shown in 
figures |S1 and IHl In these figures, the excitation proba- 
bility is plotted as a function of the change rate g and 
as a function of the particle-number N respectively. The 
dots represent the results from numerical simulations and 
the solid lines represent evaluations of the analytical for- 
mula (UI). It is seen that the numerical results for the 
excitation probability scale linearly with the number of 
particles N and they scale quadratically with the change 
rate g, as it is stated by formula {Tj). 
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VI. STABILITY 

In this section, the imphcation of perturbations that 
accompany the adiabatic process shall be discussed. 
These perturbations occur, for example, if the adiabatic 
process is simulated on a quantum computer and quan- 
tum gates are not implemented perfectly. The problem 
that arises is that perturbations usually break the sym- 
metry of the original Hamiltonian, such that, in princi- 
ple, excitations of levels with very low energy may occur. 
The question that has to be answered in this connection 
is whether these excitations spoil previous results or not. 

In the case of the quantum Ising model, the answer 
to this question can be found using time-dependent per- 
turbation theory with respect to the perturbation. The 
perturbation is thereby assumed to be given by 

If the adiabatic process consists in the slow decrease of 
the field-strength g from go ^ 1 to 51 < 1, then, after 
crossing the critical point, this perturbation leads to a 
strong excitation of the (asymptotically degenerate) first 
excited state and a very small excitation of states with 
higher energy. Nevertheless, the mean energy of the sys- 
tem after excitation is still very low compared to the 
width of the spectrum Ai?, namely, 



{H) - Eq 
AE 



const e 



This term doesn't depend on the number of particles N 
and is negligible for e <g; 1 even for large N . Thus, it 
can be concluded that, even though an adiabatic process 
accompanied by small perturbations will not leave the 
system in the ground state, it will prepare the system in a 
state with very low temperature. Such a low-temperature 
state also shows properties which are interesting to be 
investigated with a quantum computer. 

In the case of the anisotropic Heisenberg model, the 
perturbation term is assumed to be more general, namely. 



V 



N-l 

3=0 



with Ej e M and \n'^^^ \ = 1. The adiabatic process, how- 
ever, is assumed to be restricted to the regime where 
the external magnetic field is much stronger than the 
internal interactions. The probability that the perturba- 
tion V causes an excitation during this process is then 
estimated as 



Pv 



where e is defined as 



< 



e 

2^ 



N 



\ 



N-l 



,0) 



From this formula it can be read off that the excitation- 
probability Pv scales polynomially with the error- 
parameter e and the number of particles N . Thus, small 
perturbations during the adiabatic process will not have 
severe implications. 



VII. CONCLUSIONS 

Summing up, adiabatic processes have been investi- 
gated in the light of the simple, exactly solvable quan- 
tum Ising model and the more complicated anisotropic 
Heisenberg model. The adiabatic processes were assumed 
to consist in the slow variation of the strength of the mag- 
netic field that environs the spin-chains. 

In the case of the quantum Ising model, the inves- 
tigation could be performed in detail and analytic re- 
sults were obtained even for processes that cross the 
quantum critical point. In the case of the anisotropic 
Heisenberg model, adiabatic processes were studied in 
regimes amenable to perturbative treatment. In both 
cases, the duration of the adiabatic processes turned out 
to scale polynomially with the number of particles the 
spin-systems consist of. 

The results that were obtained are relevant for the 
treatment of bosons in optical lattices and for Adiabatic 
Quantum Computation. In Adiabatic Quantum Compu- 
tation, light is shed on the efRciency of adiabatic quan- 
tum algorithms that investigate the ground state of spin- 
systems. In the treatment of optical lattices, insight is 
gained into the dynamics of quantum phase transitions. 

As an outlook on future work, adiabatic processes in 
more complicated, yet unsolved models, such as higher- 
dimensional spin-systems or spin-glasses, could be stud- 
ied, such that information is gained about quantum phase 
transitions in these models and about the efRciency of 
adiabatic quantum algorithms investigating the ground 
state of these models. 



APPENDIX A: DIAGONALIZATION OF THE 
QUANTUM ISING-HAMILTONIAN 

Jordan- Wigner transformation J£] of Hamiltonian lfT|l 
leads to a Hamiltonian consisting of two parts: a part 
that is quadratic in terms of Fermi operators and a non- 
quadratic part: 

N-l 



H/J = -J2 + 44+1 + h.c.) - 2.g4c,} + 

i=0 

+ (c]v_iCo + c]v„ict + h.c.) (e*"* + l)-9N 



The Fermi-operators Cj are chosen to be periodic, i.e. 
Cj+AT = Cj. The symbol N denotes the fermion num- 
ber operator which is defined as N 



T^-\^c- 



Even 



though the non-quadratic term can be neglected for many 
calculations because it only makes changes of order to 
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eigenvalues and eigenstates it is fundamental for 

the investigation of the quantum dynamics in the vicinity 
of the critical point 5 = 1. This is because the ground 
state and the lowest excited state approach according 
to ^ in the vicinity of the critical point as the number 
of particles N increases (see figure Thus, corrections 
of order -i- cannot be neglected even if N is large. 

Since the non-quadratic term is invariant under 
linear transformations between Fermi-operators [T^ . di- 
agonalization can be performed separately in two sub- 
spaces consisting of an even and an odd number of 
fermions. In the odd-fermion-number subspace, e^'^^ 
simplifies to —1 and the remaining quadratic Hamilto- 
nian can be diagonalized by linearly transforming the set 
of Fermi-operators Cj into a set of Fermi-operators fjj in 
terms of which Hamiltonian Q is diagonal: 



H/J 



N-l 
k=0 



(Al) 



The energy A^ of a single 77fe-fermions is calculated as 
g-l, k = 0) 



Afc 




2gcos^, k^O 



(A2) 



In the even-fermion-number subspace e*'^^ equals 1 and 
diagonalization yields 



H/J = gN+J2Aj4vk-^ 



k=0 



with 



A. = 2Wl 



9 o 27rfc 
.g2 + 2g cos — 



(A3) 



(A4) 



being the energy of one single 77/j-fermion. 

The ground state of Hamiltonian ^ is the vacuum- 
state I vac) in the even-fermion-number subspace. The 
first excited state is the 1-fermion state 77^ | vac ) lying in 
the odd-particle subspace. The energy difference between 
these two states amounts to 



Ai?(5) = (.9-l) + y(.9-l)2+g(^- 

for ^ 1. This energy difference tends to zero for 
g < 1 in the thermodynamic limit, which is known as 
asymptotic degeneracy p^. 



APPENDIX B: CALCULATING THE DYNAMICS 
OF THE ISING MODEL 

Adiabatic processes are conveniently treated by means 
of the Adiabatic Approximation 23] . The Adiabatic Ap- 
proximation allows to calculate in first order the prob- 
ability that, during an adiabatic passage, states above 



the ground state are excited. If it is assumed that the 
adiabatic process consists in the variation of the parame- 
ter s from to 1 during time T, the first order excitation 
probabilities read 



Po^.(T)<^ max 



ujjo{s) 



(Bl) 



with 



^^■"^'^ ^ 

fiwjo(s) = Ej{s) - Eo{s). 

I V'j(s) ) and Ej{s) thereby denote the j*'' Eigenstate of 
H{s) and its corresponding Eigenvalue. A definite state- 
ment about the duration of the adiabatic process is finally 
obtained by demanding that the overall excitation prob- 
ability must be negligible, i.e. pe = J2j^oPo->jiT) < 1- 
In the adiabatic process investigated here, the field- 
strength g is varied linearly from go to gi during time T. 
In terms of the previously introduced formalism, g is de- 
pendent on the parameter s according to the formula 

g{s) = go + s{gi - go) 

and s is varied from to 1 as time evolves from to T. 

Since the adiabatic passage starts from the ground 
state (lying in the even-fermion-number subspace) and 
the evenness and oddness of states is conserved during 
time evolution because of the Z2-symmetry, calculations 
can be restricted to the even-fermion-number subspace. 
In this subspace, the dynamics are governed by Hamilto- 
nian ((X3|) . 

The excitation probability can now be calculated by 
means of formula However, the determination of 

the functions Ajo(s) requires the knowledge of the deriva- 
tive of H{s) with respect to the parameter s. This deriva- 
tive can be calculated easily if it is taken into account 
that (for an odd number N of spins) 



ds 



with 



sm 



27rfc 
N 



dg 

rfITT5(s)2 + 2g(s), 



2irk 
N 



The derivative H'(s) couples the vacuum state | vac{s) ) 
only to the vacuum state itself and to the two-fermion 

states 



\nis)) :^r]l{s)i]ljs)\vac{s)), n =!,...,[ 



N-1 



-J- 



The only states that are excited (in first order) during 
the adiabatic passage are therefore the two-fermion states 
I n{s)). The matrix element of H'{s) that describes the 
coupling between the states | n{s) ) and | vac{s) ) reads 



[n{s) \H'{s)\vac{s) 
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In addition, the energy difference between | n(s) ) and the 
ground state amounts to 

hujnois) = 2JA„(s). 

The probabiUty that the two-fermion states | n{s) ) are 
excited can now be calculated using formula (jBljl : 



Po- 



< 



1 



max 



24 J2 T2 0<s<l 



Us) 



A„(s) 



The overall excitation probability p^;, i.e. the probability 
that any state above the ground state is excited during 
the adiabatic process, is equal to the sum over all terms 
Po^n- This sum can be evaluated for N ^ 1 either by 
replacing the sum by an integral yielding result Q or by 
explicitly summing over approximated addends yielding 
results © and 



APPENDIX C: CALCULATING THE DYNAMICS 
OF THE HEISENBERG MODEL 

As in the previous section, the probability that the sys- 
tem is, at the end of the adiabatic process, in an excited 
state can be calculation by means of formula (|Bip . The 
calculation of the quantities Ajo{s) and fkOjols), how- 
ever, requires the knowledge of the eigenstates and eigen- 
values of Hamiltonian The series expansions about 
5 — *■ oo of the ground state energy Eo{s) and the ground 
state I ipo{s) ) read 



Eojs) 
JN 



and 



-gis) + 



|Vo(s))(x|^(°)> + -i^|^(^^)+0 



1 2 



The coefficients | (/?q°^ ) and | ipg''' ) are thereby defined as 



(1) 



I (0) ' 

I (1) ' 



N~l 



i(A,-A,)^|j,j + l 



J=o 



The symmetries that the ground state possesses are the 
Z2-symmetry 1 and the translational-symmetry fc = 0. 
The only states that can be excited during a time evolu- 
tion according to H(s) are states that possess the same 
symmetries as the ground state. The lowest states that 
meet this requirement are the states with group-index 
n = 2 and translational-symmetry k = 0, namely, 

\i'20d{s)). 

The energies of these states are £'20d(s). In the case of an 
odd number of particles TV, the series expansions about 



.g — > oo of E20d{s) and | V'20d(s) ) read 



E2ad{s) 
JN 



-5(.)(1-^ -A. 1-- 



and 



\ip2Qdis)) (X l^fad) 



Aylcos — 



27rd 



O 



'P20d / 



o 



rin2 



The index d may thereby assume the values 1 , . . . , ^ ^ ^ ■ 
The coefficients | ^^20d ) projection of | v^joi ) on 



are given by 



I (0) > 



I (1) > 



N-1 
2 



1 



= / sm 
1 ^ N- 



2ndj 



- 20j , 



(A. 



N 



2-Kd 



N + 1 N + 1 



The kets | T20J ) that appear in the upper formula denote 
the 2-particle states with translational symmetry fc = 0. 
They are defined as 



-20j 



i + d). 



In order to calculate the functions Ajo(s), the matrix 
elements of the derivative H'{s) between the ground state 
I V'o(s) ) and the excited states | ip20dis) ) must be deter- 
mined. It is advantageous to abbreviate these matrix 
elements by the notation Ti.da{s): 

Hdois) {,P20d{s)\H'{s)\Ms)) 
The derivative H'{s) is equal to the expression 



H'is) = -Jg'{s) 



N-l 
i=0 



Using this expression and the series expansions of | iPq{s) ) 
and I V'20(i(s) ), it is straightforward to determine the ma- 
trix elements Ti.do{s)- The contribution to Tido{s) to O*''- 
order in - is zero. This is because the derivative H'(s) 

does not couple the 0*'*-order-coefficients | ip^^^ ) and 
I ^^^l^). To order in \, the terms (<^« \E\s)\ ) 

and (</'2o!i l-^'l*)! ¥'0^'' ) contribute. Higher-order contri- 
butions to Tido will not be taken into account. In sum- 
mary, the series expansion to lowest order in i of the 
matrix element Ti^o reads 

^'^"(■'') = ^ { (^2oil^'(^)l4"^) + 
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A simplification of this expression yields 

. X , ^ / N . 2nd ^ ri" 

T^^o s = -2(A^ - A„) ^^H^W s n hO - 

''"^ ' ^ " yg{s)\lN+l N+1 [g 

The probabihty po^d that the states | ip20d{s) ) are ex- 
cited can now easily be calculated using formula (|Bip . 



The overall excitation probability p^; is equal to the sum 
over all terms po^d- In the limit of a large number of 
particles {N ^ 1), this sum can be approximated by 
an integral. This approximation leads to a simplified ex- 
pression for the overall excitation probability , namely, 
result (71 . 
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